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Abstract 

In this manuscript we present a detailed investigation of the form factors of bound- 
ary fields of the sinh-Gordon model with a particular type of Dirichlet boundary 
condition, corresponding to zero value of the sinh-Gordon field at the boundary, at 
the self-dual point. We follow for this the boundary form factor program recently 
proposed by Z. Bajnok, L. Palla and G. Takacs pQ, extending the analysis of the 
boundary sinh-Gordon model initiated there. The main result of the paper is a con- 
jecture for the structure of all n-particle form factors of the boundary operators d x <j) 
and (d x (j)) 2 in terms of elementary symmetric polynomials in certain functions of the 
rapidity variables. In addition, form factors of boundary "descendant" fields have 
been constructed. 
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1 Introduction 



In the context of 1+1-integrable quantum field theories (QFTs) form factors are defined as 
tensor valued functions, representing matrix elements of some local operator 0{x) located at 
the origin between a multi-particle in-state and the vacuum 

F^-^(e 1 , ...,e n )-.= (010(0)^(^)^(02) . . . z, n {e n )\Q). m . (1.1) 

Here |0) represents the vacuum state and Z^.(9i) are vertex operators which provide a general- 
ization of Bosonic or Fermionic algebras and allow for the definition of a space of physical states. 
They carry indices /ij, which are quantum numbers characterizing the various particle species, 
and depend on the parameters 9i, which are called rapidities. In 1+1 dimensions the energy and 
momentum of the particles can be expressed in terms of the rapidity variables as 

Ei = rriiCOshOi and Pi = mjsinh#j, (1-2) 

where rrii is the mass of the particle considered. The braiding relations amongst the operators 
Z^.(9i) and their Hermitian conjugates are known as Faddeev-Zamolodchikov algebra For 
theories without back-scattering (diagonal theories) they are given by 

Zfiii&^Z^^j) = S^^Oi^Z^^B^Z^Oi), (1.3) 

zlmz^e,) = s^e^zl^zlm, (1.4) 
z Mi (^)4,(^) = Sw(-0 ij )Zl.(0 j )Z lk (e i )+2i r 6 IHflj Wi)> ( L5 ) 

where S^ jfJij (8ij) is the 2-particle scattering matrix and we introduced the variable Oij = Q{ — 8j. 

The form factor program for integrable models was pioneered by P. Weisz and M. Karowski 
[3J0I] and thereafter developed to a large extent by F. A. Smirnov, who also formulated some of 
the consistency equations for form factors [§]. In these seminal works, the fundamental properties 
of form factors in 1+1-dimensional theories were established (see [Z] for a recent review). It was 
found that the form factors of local operators can be obtained as the solutions to a set of 
consistency equations whose origin is based on physically-motivated requirements. The solution 
to these equations allows in principle the computation of all n-particle form factors associated 
to any local field of the massive QFT. Following this program, form factors of a large class of 
integrable models have been computed (see e.g. |l]-|29j and references therein). Once the form 
factors of a certain operator are known they can be used for many interesting applications like 
the computation of correlation functions (see e.g. [S], and HUIISl)) the identification 

of the operator content of the perturbed conformal field theory (CFT) (see e.g. EDI EI] and 
HH-EZ1) an d the explicit computation of quantities which characterize the underlying CFT (see 
e.g. [S21 OH EZl EH and H3-I22])- 

So far, we have briefly recalled the program for bulk integrable QFTs, that is QFTs without 
boundaries or defects. Given that most realistic physical systems possess boundaries and/or 
defects, it is both interesting and natural to ask the question whether or not the integrability 
of certain 1+1 dimensional QFTs may be preserved in the presence of boundaries or defects. 
The study of the conditions under which that question can be answered in the affirmative has 
attracted a lot of attention in the last two decades, with many works devoted to the study of 
integrable QFTs with boundaries [S3"|SH1 and defects (1231112 1^- A natural further step is 
trying to extend the bulk form factor program just recalled to boundary and defect theories. 
Here we want to concentrate on the boundary case. Given an integrable QFT with a boundary 
there are two possible approaches to the problem of computing correlation functions: 
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Boundary at t = 0: we may place the boundary at t = extending along the space direction 
(t represents the time coordinate) . In this case it is natural to characterize the presence of the 
boundary by means of a boundary state \B). The Hilbert space is exactly the same as for the 
bulk theory, and therefore form factors can be completely characterized by those in the bulk. 
The only additional information needed is an explicit description of the boundary state, which 
has been provided by S. Ghoshal and A.B. Zamolodchikov 36J, 

\B) :=exp(^ R ltl (^-B)Z fH (-0)Z IH (9)d6\ |0>. (1.6) 

Here R^. (9) is the reflection amplitude off the boundary, which we take to be diagonal in this 
case (particles just reflect off the boundary with a certain probability, without changing their 
species). Expanding the exponential above, boundary form factors can be expressed in terms 
of bulk form factors, although the former will then be an infinite sum of the latter. This sum 
is in general very involved, but can be performed analytically for some correlation functions of 
free theories |42] . This picture can be generalized to the defect case, by replacing the boundary 
state above by a defect state. A realization of this state and explicit form factor computations 
have been carried out for free theories |391 14U| with different types of defects. More recently it 
has been proven in that a defect theory may always be regarded alternatively as a boundary 
theory by means of the so-called "folding trick". Exploiting this correspondence, the authors 
also constructed the defect state for integrable models with purely transmitting defects. 

Boundary at x = 0: we may place the boundary at x = extending in the time direction (x 
represents the position coordinate). In this case the Hilbert space will change with respect to 
the bulk only half of the physical states remain linearly independent in the presence of 

a boundary. This is so since in addition to (jl.3)) - (jl.5j) the further constraint 

Z^) = R^Z^i), (1.7) 

holds. Recently, a generalization of the bulk form factor program for theories with boundaries 
and for boundary operators (that is, operators defined on the boundary) has been proposed 
PQ, and its validity tested against previously known results for various QFTs. It has been 
established that the relations ([1.3)1 - (|1.5|) and ([1.7)1 together with the analytical properties of S- 
and i?-matrices lead to a set of consistency equations for the boundary form factors, which can 
be solved for particular models along similar lines as for bulk theories. Some of these solutions 
were found in pQ, although in most cases only for low particle numbers. The main aim of this 
paper is to extend some of the results in 1 by carrying out a detailed study of the boundary 
form factor solutions of the sinh-Gordon model with a particular type of Dirichlet boundary 
conditions (corresponding to vanishing value of the sinh-Gordon field at the boundary) at the 
self-dual point. 

The paper is organized as follows: In section 1 we review the boundary form factor program 
as proposed in pQ. In section 2 we recall the S- and /^-matrices of the sinh-Gordon model with 
Dirichlet boundary conditions and vanishing boundary field. We also present the form factor 
equations found in for this particular model. In section 3 we re-write these equations in terms 
of elementary symmetric polynomials and analyze their special properties at the self-dual point. 
In section 4 we provide new explicit solutions to the form factor equations up to 16-particles 
and formulate a conjecture for the form of all n-particle form factors of the boundary operators 
d x 4> and (d x 4>) 2 . In section 5 we find solutions to the form factor consistency equations for other 
boundary fields. We summarize our conclusions in section 6. Finally, we list various explicit 
form factor formulae in the appendix. 
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2 The boundary form factor program 

In this section we will review the main aspects of the boundary form factor program proposed in 
for diagonal scattering theories possessing a single particle spectrum and no particle bound 
states (such as the sinh-Gordon model). For such a model, the particle indices fa in IJl.lj) can 
be dropped and the boundary form factor axioms written as 

F n . . . , 9i,0i + i, . . . n ) = S(8i i+i)F n (0\, . . . , 0i + \, 9i, . . . n ), (2.1) 

-^n (01> • • • 3 ^n-l) On) = Rn(8n)Fn (0\, . . . ,0 n -i, —0 n ), (2-2) 

F°(9 u 9 2 ,...,9 n ) = R n (in-e 1 )F°(2m-6 1 ,e 2 ,...,0 n ), (2.3) 

The first two axioms follow from the braiding relations ()1.3|) and Q1.7JI . whereas the last axiom 
expresses crossing symmetry in the presence of a boundary. In addition we have the kinematical 
residue equations which generalize the corresponding ones in the bulk 

- iRes- eo=eo F° +2 (0 o + in, Mi • • • , n ) = \1 - f[ S(0 oi )S(0oi)\ F%{0 X , . . . , n ), (2.4) 

with 0ij := 0i + 0j. It was argued in 1. that a second kinematical relation which has no 
counterpart in the bulk form factor program emerges whenever the reflection amplitude R(0) 
has a pole at = wr/2. We will not report this equation here, since it is trivial for the particular 
model we want to study later. 

As in the bulk case, further recursive relations exist in the presence of particle bound states 
but we will not review them here. Mimicking the solution procedure developed for the bulk 
case and diagonal scattering matrices, an ansatz for the solutions of the above equations can be 
made, which solves axioms (|2,l|) - (l2.3j) by construction and at the same time makes explicit the 
form factor pole structure 

F^0 1 ,...,0 n )=H^(y 1 ,...,y n )f[r(0 i ) ft IMMlA. (2 . 5 ) 

i=l l<i<j<n W 

Here yi = 2cosh6>j, H® is a constant and Q^ivi, ■ ■ ■ ,Vn) is an entire function of the variables 
yi,...,y n . The function f(0) is the same minimal (2-particle) form factor found in the bulk 
case and r{0) is the minimal 1-particle form factor. Integral representations for these minimal 
form factors can be readily obtained from the representations of the S'-matrix and reflection 
amplitude R(0), respectively. 

3 Boundary form factors of the sinh-Gordon model with Dirich- 
let boundary conditions 

The particular example we want to study in this letter is the sinh-Gordon theory with a particular 
type of Dirichlet boundary conditions. The sinh-Gordon model is one of the simplest integrablc 
models possessing a non-trivial scattering matrix since its spectrum consists of a single particle 
and no particle bound states occur. The bulk scattering matrix |441 1451 HH] is a pure CDD factor 
given by 

S(0) = -(-B)e(B-2) e , (3.1) 



in terms of the blocks 



[X)d ~ sinh \{0-*f)' { } 
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The parameter B £ [0,2] is the effective coupling constant which is related to the coupling 
constant (5 in the sinh-Gordon Lagrangian |47| I48j 



1 ,„ ,,o m 
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C=^-(d^f -^cosh(/3c/>), (3.3) 
1 P 



where m is a mass scale and (f> is the sinh-Gordon field, as 

2j? 
8vr + /3 : 



under CFT normalization 3 . The S-matrix is obviously invariant under the transformation 
B — > 2—B, a symmetry which is also referred to as week-strong coupling duality, as it corresponds 
to B(P) — > B((3~ 1 ) in (|3.4fl . The point B = 1 is known as the self-dual point. When an integrable 
boundary is introduced and Dirichlet boundary conditions imposed the corresponding reflection 
amplitude R{9) was found to be [IH] 

This expression was obtained by analytical continuation of the sine-Gordon breather reflection 
amplitudes calculated in [SO]. In (|3.5f) E is a continuous parameter which is proportional to the 
boundary value of the sinh-Gordon field as 

p 

At the special value E = 0(0) = which we will consider in this manuscript, the reflection 
amplitude (|3.5|) has no pole at 9 = in/2 and therefore the ansatz (|2.5|) captures the full pole 
structure of the boundary form factors. In addition, since Dirichlet boundary conditions imply 
a fixed value of c/>(0), the boundary operator content of the theory consists solely of products 
of space-derivatives of the boundary field like d x cp, (d x (p) 2 , that is, derivatives in the direction 
perpendicular to the boundary. The recursive equation obtained by substituting (|2.5|) in the 
kinematic residue equation ()2.4|) was derived in and takes the form 

Qn+2(-y>y>yi> ■■■,y n ) = -P n (y,yi, ■ ■ ■ ,yn)Q°(yi, • • ■ ,y n )- (3.7) 

The polynomial P n was also given in pQ. Here we present it in a slightly different form, in terms 
of elementary symmetric polynomials in the variables y±, . . . , y n , 

P n (y,yi,...,y n ) = l ^ U+U} -> £ (-^_)- fc (-^ + )"f a^sin^ -p)/2), (3.8) 
^ k,p=0 

which are defined through the generating function* 

n n 

Y[{ z + yk ) = Y J Z n - k ot, (3.9) 



k=l k=0 



The variables u± are given by 



u± = 2cosh(6 ±iirB/2). (3.10) 



*In what follows we will employ the simpler notation <r£ = a^, since the number of variables n on which the 
polynomials depend is usually obvious from the context in which they occur. 
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3.1 Recursive relations at the self-dual point 



In this section we initiate the presentation and discussion of the new results obtained in this 
manuscript. Whereas in computations up to n = 5 particle form factors and for generic values 
of B where carried out, here we have concentrated our analysis on the special case B = 1, that 
is the self-dual point. In this section we will show how at this point a new factorized structure 
is found, both of the form factor equations and their solutions, which will allow us to find form 
factor solutions up to much higher particle numbers (n = 16) and to put forward a plausible 
hypothesis for the structure of the solutions of the form factor equations for arbitrary particle 
numbers. Let us start by noticing that for B = 1 the variables u± become proportional to each 
other 

u± = ±2isinh0, (3.11) 
and this allows us after some simple algebra to bring P n into the factorized form 

P n (y,yi, . . . ,y„) = P%(y,yi,.- ■ ,y n )P£(y,vi, ■ ■ -,y n ), (3.12) 

with 

n 

PZ(.V,Vi,...,y n ) = (-iu + )*Wj2o%(-iu,+)-r cos(7rp/2), 

p=0 

n 

p°( y , yi ,...,y n ) = (-ir + \-iu + fi^}- i j2*pH"+r p MKpm, (3.13) 

p=0 

where [x] denotes the integer part of x. Clearly, the polynomials P% and P° involve only elemen- 
tary symmetric polynomials of even and odd degree, respectively. In addition, the coefficients 
of the symmetric polynomials in P% and P° are always even powers of — iuj+, that is, P% and P° 
are polynomials in (iuj + ) 2 = y 2 — 4, of which the first few examples are 



PI = 1, Pf = ax, 

Pi = (iu; + ) 2 -a 2 , P 2 ° = -a u 

Pf = {iuj + ) 2 -a 2 , P 3 ° = (iuj + ) 2 ai-cr 3 , 

PI = (iw+) 4 - (iu + ) 2 a 2 + 0-4, P 4 ° = -{iu + ) 2 <Ti + a 3 . 



(3.14) 



It was noticed in that for E = the symmetry of the Lagrangian under <j) — * —(f), where <f> 
is the sinh-Gordon fundamental field, is preserved in the boundary theory. This means that, 
depending on the operator O considered, only the even or odd particle form factors will be 
non- vanishing, exactly as in the bulk theory [131 1141 129j . 

In the first non-vanishing form factors of two operators identified as d x (p and (d x (j)) 2 were 
computed. Setting B = 1, the solutions found become 

Ql x4> = 1, (3.15) 

Q (d x <pT = (316) 

Qf* = -<n, (3-17) 
Qf^) 2 = a 2 (4 + a 2 ), (3.18) 
Qfr* = fJl (4 f T 1 +a 3 )(4 + a 2 ). (3.19) 

Here we have presented the solutions for both operators in order of increasing particle number 
to make certain features more apparent: first of all, all Q n factorize into a factor containing only 
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elementary symmetric polynomials of even degree and another factor containing only those of 
odd degree. Secondly, some of these factors repeat in various form factors (like 4 + a<i which 



appears both in Q^f x ^ and Q^ 9 ), so that some sort of pattern seems to emerge. In this 
letter strong evidence will be provided that these two features are actually general: all functions 
Qn*, Qri x ^ factorize and a general pattern for the n-particle solutions exists, such that once the 
form factors of d x (j) are known also those of (d x (j)) 2 are automatically determined and viceversa. 

Recalling that the P-function (|3.8j) factorizes in a similar way, our results suggest that 
the recursive equations (|3.7j) can be split up into two simpler recursive equations involving 
elementary symmetric polynomials of even and odd degree, respectively 

Q°( yi , ...,y n )= Q° le (yi, y n )Q° l °(yi, • • • , y n ), 



with 



-y,y,yi, 
-y,y,yi, 



,y n ) 

,Vn) 



P5(v,vi,-- 
-PZiv,vi, 



,y n )Q° le (yi,.. 
■ ,yn)Q° l °(yu 



,Vn) 



,y n ) 



(3.20) 



(3.21) 



In this new notation, the solutions above become: 



eft* 

(d x <f>) 2 



Q 

^3 
^5 



1, 

-1, 
1, 

4 + CT 2 , 
-(4 + (T 2 ) 



eft* 

^3 

Qd x <j> 
5 



1, 



-o-i, 



(3.22) 



(7 



1; 



-01(4(7! + (T 3 ). 



This factorization simplifies greatly the process of finding solutions of Q3.7|) for higher particle 
numbers. From (j3.21j) and (|3.13j) two important characteristics of the solutions of (|3.21|) can be 
deduced: 

• the functions Q^f^ an d Q^nfi ^ always coincide, but for the number of variables entering 
the symmetric polynomials and a sign. The reason is that PJ^ an d fJn+l are m ^ ac ^ ^ ne 
same in terms of elementary symmetric polynomials and that the input values Q^ x *^ 6 = 
—Q^ e = 1 are the same but for a sign. 

• Similarly, Q^n+l '° an d Qtn+2 a ^ so coincide but for an extra —u\ factor in Q2n+2- This is 
again due to the fact that the corresponding P-functions, P% n+ i and are identical 
up to a sign in terms of elementary symmetric polynomials and that the initial values 



2 ~ 

631). 



<j\ and Q\ 



(a^) 2 |o 



1 differ only by a u\ factor which satisfies the first equation 



m 



These features can be found on the few examples listed above. Most importantly, they allow 
us to conclude that, at the self-dual point, the knowledge of all non- vanishing form factors of 
the operator {d x (f)) 2 automatically implies the knowledge of all non-vanishing form factors of the 
operator d x cp through the relations 



n {dx$) 2 \e 
^2n 

n (d x <t>) 2 \o 
V2n 



-Q 



d x <p\e 
2n+l> 



-<?lQ 2n -l 



for n € Z"* 
for n € 



(3.23) 
(3.24) 



Note that these equalities are only true when the functions are expressed in terms of elementary 
symmetric polynomials. The latter will obviously depend on different particle numbers on the 
r.h.s. and l.h.s. of the equations (and therefore, on different numbers of variables). 
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4 Solutions for higher particle numbers 

We shall now attempt to find new solutions to (|3.21f) in terms of elementary symmetric polynomi- 
als in the standard way employed previously for several bulk theories (see e.g. [HAHHHm 1221125] ). 
This solution procedure can be rather tedious as the number of terms on the r.h.s. of IJM.21|) 
increases rapidly with the particle number. Prom these terms one must be able to "reconstruct" 
the l.h.s. of the equation when the two first variables are not y and — y but generic. This recon- 
struction can be performed by employing the following relations between elementary symmetric 
polynomials depending on n + 2 and n variables 

<rk(-y,y,yi,---,y n ) = -y 2 ^k~2(yi,---,y n ) + a k {yi,... ,y n ), (4.1) 

with a -i = and ctq = 1. The relations (|4,lj) allow us to make a guess about the possible 
origin of each of the terms on the r.h.s. of 1)3.21)1 . For example, if a term appears it can 
only originate from a term a^ +2 in the l.h.s. and if a term — y 2 cx™ appears, then it can originate 
either from a ^+2 or a a 2 +2<T n + ' 2 term on the r.h.s. By consistently accounting for each term 
on the r.h.s. of ()3.21j) a unique polynomial solution to the equations can be found. Solving the 
recursive equations systematically in this way we have found explicit solutions up to 16-particle 
form factors. All the solutions found exhibit the following general structure: 

Q2n - y~ i ) V 2 n-2 /^n-2 + °2n-«2nJ l^.I) 

n {d x <j>) 2 \o _ i -,\n n {d x 4>) 2 \o j_JH„ po (a o\ 

Q2n - ^2n-2 ^2n-3 + ^1^2n-l-«2n> ( 4 "JJ 

Let us now explain the various terms appearing in 1)4.2(1 - 1)4.3(1 : 

• The first term on the r.h.s. of the equations shows a recursive relation between 2n- and 

(2n — 2)-particle solutions of (|3.21() . However, the polynomials Qfn-2 ^ an d Qzn-2 
appearing on the r.h.s. of the above equations are not any more the solutions of ((3.21)1 . 
They are identical to them in terms of elementary symmetric polynomials, but these are 
now polynomials on all 2n variables involved on the l.h.s. of the equation. 

• The functions fa are defined as follows: 

fi = 1, fj, 1 = a 1 , 

Hi = 4/^ + (7i, for i>l, (4.4) 
in terms of elementary symmetric polynomials. Equivalently 

[i/2] 

fJ, i = J2 4k<T i-2k, (4.5) 

where, as before [x] stands for the integer part of x. The emergence of such particular 
combinations of elementary symmetric polynomials is rather natural given their reduction 
properties, which can be deduced from ((4.1)1 : 

Vk{-y>y>yu--->yn) = (4-y 2 W_ 2 (yi>--->yn) + ^k{yi,---,y n ), (4.6) 

where fj,_i = and n = 1. As we can see the variable y appears always in the combination 
4 — y 2 , which is precisely the same variable entering the polynomials P e and P°. 
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• Further, the conjectured expressions (|4.'2j) - (|4.3|) involve the additional terms a2nR% n an d 
"i cr 2n-i^2n' respectively. One can easily convince oneself from the previous definitions, 
that these are in fact the only terms containing the elementary symmetric polynomials a 2n 
and a 2n -i, respectively Therefore a 2n R\ n anci a \ (y 2n-\R 2n are the terms containing the 
elementary symmetric polynomial of maximum degree. This definition makes it straight- 
forward to reconstruct the functions i?| n and K^n f° r particular cases up to high particle 
numbers. From (|4~2|) . (J4~3|) . (jS2U>, ff^TTf) and it follows that they must satisfy the 
following reduction properties: 

-y 2 Rtn{-y, y, 1/1, ■ ■ ■ , yzn-zj^zl = PL-zQfn-l |e (yi> ■ ■ ■ > 2/2^-2) 

+(-!)" Q&_f |e (-y, y, Ift, . . . , 2/2n- 2 )((4 - y 2 )^l4 + (4-7) 

y\al n - 2 fR° 2n (-y,y,yU ■ ■ ,&- 2 )^3 = P2n-2Q { 2nT°(yi, ■ ■ -,V2n-2) 

+(-l) n Q&_f |o (-y, y, yi, ... , Mta- 2 )((4 - y 2 )^: 2 + a 2 ™: 2 ), (4.8) 

where, in order to avoid confusion, we have made the variable dependence explicit and 

a 2 k n ~ 2 = a k (y t,..., y 2n -2) and /u 2n ~ 2 = n k {yi, . . . , y-in-i), (4.9) 

for any values of k. Notice that, once the conjectured formulae (|4.2|) - (|4.3[) are assumed 
to hold, i?| n and R\ n are the only unknowns, which means an enormous simplification of 
the original problem. Explicit formulae for R\ n and R^n up to n = 8 can be found in the 
appendix. 

(d <t>) 2 \e 

In addition to the structure just outlined, a simple recipe allows to relate the solutions Q 2n 
andQg^ 2 ' . Thelatt er can always be obtained from the former by performing the replacements 
Cfc — > ffc+i and \x k —>■ fJ, k +i and multiplying by a global factor —a\. Once these transformations 
have been done, we must still introduce as many factors cx\ as to achieve that each term in 
Q2n+2 '° i s exactly a product of n + 1 symmetric polynomials. For example, if a term erg would 
appear in Q^ x ^ e , then it would imply a term —<j\g\g§ in 

g(jW) a |e_ ThuS; provided 

we know 

Q 2 n ' 6 we can systematically determine Q 2 j^ '°- From these solutions, as concluded above 
also all functions Q^n+i and Qtn+l can be obtained. 



5 Solutions for other boundary fields 

In this section we will address the problem of finding solutions of the form factor equations 
for fields other than d x (j) and (d x (j)) 2 . It is worth recalling that the fields d x cj) and (d x cj)) 2 were 
identified in 1 by analyzing the asymptotic behaviour of their form factors. Although this 
identification seems plausible, it would be very interesting to strengthen it with some numerical 
work. 

A particular way of looking for new solutions of the form factor equations is to find entire 
functions • • • , y n ) such that if F®(yi, . . . , y n ) is a solution, then the product 

F°\yu ...,y n ) = r n ( yi ,..., yn)F°(y u ...,y n ), (5.1) 

is also one. In the bulk case, such type of form factors (with variables yi = e 8i + e~ Si replaced 
by X{ = e Si ) have been related to "descendent" fields for the first time in the work of J. L. Cardy 
and G. Mussardo |5j for the Ising model. The name descendent must be understood in the sense 
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that if O is related to a spinless primary field of the underlying CFT in the ultraviolet limit, 
then O 1 would be related to a descendent field of O in the same limit. Which particular class of 
descendent fields O' relates to is determined by the scaling properties of the function 1^ under a 
rapidity shift, that is, the spin of O' which we call s. Solutions to the bulk form factor equations 
of the type ()5.1|) have been also found for the sinh-Gordon ^3] an d Yang-Lee [SJ models. 

In our case, the link between form factors of the type 1)5. 1J) and descendent fields in the 
conformal limit is a priory less clear, as the operator content of the CFT at the boundary is 
not classified in terms of spin representations. However, we may still attempt to find solutions 
of the form 1)5.1)) . They will be related to new fields of the QFT, that is higher derivatives 
of <j) which is very natural to relate to descendent fields at conformal level. In this case the 
superscript s in 1^ is just the leading order behaviour of when all rapidities tend to infinity. 
The analysis is in fact completely analog to the one presented in ^3] for the bulk sinh-Gordon 
model. We just need to replace the variables Xi by the variables j/j. First of all, it is trivial 
to see that the form factor 1)5.1)1 will automatically satisfy all consistency equations (|2.1[) - (|2.3|) 
provided I^(yi, ■ ■ ■ > Un) can De entirely expressed in terms of elementary symmetric polynomials 
in the variables {yi, . . . , y n }. Secondly, if ()5.1)) is to satisfy the kinematic residue equation 1)2.4)1 
then 1^ must satisfy the following reduction property: 

!n+2(-y,y,yh---,yn) =i£(yi,...,y n )- (5-2) 

This is exactly the same equation found in [HI EH • A basis of the space of solutions to 1)5.2)1 has 
been found to be 

Il s ~ l = (-l) s+1 detX, (5.3) 

with I being a matrix of entries 

lij = &2j-i and lij = a 2 j-2i+2, (5.4) 

for j = 1, . . . , s and i = 2, . . . , s. As in the bulk case det I is of order 2s — 1 in terms of elementary 
symmetric polynomials. 

To finish this section it is worth mentioning that it would be very interesting to establish 
rigorously a correspondence between the number of solutions to the form factor equations for 
fixed s and the number of descendent fields of the underlying CFT at level s. For this one should 
proceed along the lines of the analysis carried out in jH) for the bulk Ising model. It is however 
to be expected that such analysis becomes more involved for the present model. 

6 Conclusions and outlook 

In this paper we have studied the form factors of two boundary operators of the sinh-Gordon 
model with a particular type of Dirichlet boundary conditions, corresponding to vanishing 
boundary value of the sinh-Gordon field, at the self-dual point. We have found that both the 
form factor recursive equations ()3,7)) and their solutions exhibit distinct properties at B = 1: 
they both factorize into one factor containing only elementary symmetric polynomials of even 
degree and another term containing only those of odd degree. This factorization has allowed 
us to obtain explicit solutions of the equations up to remarkably high particle numbers and to 
provide a conjecture for the structure of the form factor solutions for arbitrary particle num- 
bers. Unfortunately we have not yet been able to provide a mathematically rigorous proof of 
that structure, although the many examples computed strongly support it. 

The finding of closed formulae for all n-particle form factors of local fields is a highly com- 
plicated task which has only been successfully carried out for few bulk scattering theories and 
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operators El QUI El El EH E3] Even when closed formulae have been found, rigorous 
proofs of their general validity have only been given in some cases |13 ( 122 ( 123]. A crucial element 
of these proofs has been the fact that the form factor solutions were given in terms of determi- 
nants whose entries were elementary symmetric polynomials. It is this determinant structure 
and the possibility of employing general properties of determinants which made these proofs 
possible. The formulae conjectured in the present paper (|4.2|l - (|4.3jl are however not of determi- 
nant form, so that general proofs should go along very different lines. In addition, it is worth 
noticing an important difference with respect to the bulk form factor program with regard to 
the structure of the form factor solutions: boundary fields do not possess a definite spin (since 
translation invariance is broken) and therefore the transformation under a rapidity shift of the 
allowed terms in the form factor formulae is not constrained. This means that not all terms in 
the functions Q n need to have the same degree. As a result, formulae tend to be more involved 
and the form factor's general structure is harder to identify. It would be very interesting to 
investigate whether or not some structures of the boundary form factors may be universal or at 
least common to many models. By this we mean structures such as the determinant formulae 
mentioned above for bulk theories. For the case at hand, it would be of course desirable to find 
a general pattern for the functions R^, i?° , which may also facilitate proofs. 

In addition, new solutions to the form factor consistency equations related to local operators 
other than d x <j> and (d x (p) 2 have been found. We have seen that the form factors of these fields 
are related to the ones of d x (j> and (d x (p) 2 by a multiplicative factor (|5.1|) . The same kind of 
structure has been previously identified for the form factors of the bulk Ising jS], sinh-Gordon ^3] 
and Yang-Lee models [H^ and has been related to descendant fields at the conformal level. It is 
natural to assume that this interpretation also holds here, however two interesting open problems 
remain: first, a numerical study of the short-distance behaviour of correlation functions of these 
fields should be carried out in order to make the identification with descendents at conformal 
level more clear. Second, a rigorous study of the correspondence between form factor solutions 
with a particular asymptotic behaviour and descendent fields for fixed level should be carried 
out along the lines of jSJ. 

We would like to conclude by saying that the program proposed in is very recent and has 
so far only been pursued for the few models treated in the original paper and for the particular 
case considered here in more detail. A lot of work remains to be done to develop the boundary 
form factor program to a degree comparable to that achieved for bulk theories. Obviously 
computing boundary form factors of other integrable models would be interesting, particularly 
so for non-diagonal theories such as the sine-Gordon model. For the latter theory boundary 
form factors have in fact been computed in [£2] by a different approach but it would still be 
interesting to employ the method of for the same computation as a consistency check and/or 
as a possible way to obtain alternative (maybe more explicit) representations of the form factors. 
It would also be very interesting to exploit these form factors solutions in similar ways as has 
been done for bulk theories: computing correlation functions, investigating the operator content 
and computing characteristic quantities of the underlying boundary CFT. Some effort in this 
direction has been made in pQ. 
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A Formulae for the polynomials R\ n and R° 2n 

For the solutions already found in and listed in (|3.22|) . it is easy to see that R\ = i?§ = 
R\ = R\ = 0. However, for higher particle form factors non-trivial i?-functions are found, 
which become more involved as the particle number increases. Expressions become simpler 
when written in terms of the polynomials: 

[i/2] [i/2] 

a* = J> + 1)4^2*, ft = -J3(fc + l)(fc + 2)4*^-2* (A.l) 

k=0 k=0 



We have constructed the following solutions^: 

R t = -1) 
R% = a 2 2 , 

R e l0 = (8(2 + a 2 )-a 4 )a fi + p 2 (al-2a 8 )+a w , (A.2) 
R i2 = -P2 a e(P<i a 6 ~ (( 16 + «2)«8 - 3ai )) + a\ - /?|(8(2/? 2 + a 4 )a 8 - 2(3 A a w + /? 2 ai 2 ), 
R\a = (8 + a 2 ) 2 a| + ag (ag - a 2 «6/34 - 4/? 2 /9i) + 2«6 (4"io - 012) /3 4 

SaiaPl ((2 + a 2 ) a 6 + W 2 Pi) + (64 (28 + a 2 (6 + /? 2 )) + af) a 8 ai /? 2 

+(4 (128 + a 2 (48 + 5a 2 )) a 4 - 2 (8 + a 2 ) - a 6 (4a 2 + /3 4 ))asaio 

+ (1792 + 4a 2 (16 (6 + (5 2 ) + a 4 ) - 12a a 6 + as) a?o + ( 8 «i2 - a i4 ) /3 2 /3 4 

+ (2a 2 ai2 + "14) «6/5 4 + 010012 ((4 - 02) «4 + 3a 6 - 4 (7a 2 + 84/3 2 - 48)) 

+ (8 + a 2 ) (a 6 - (8 + 3/? 2 ) /? 4 ) 03012 + (208 + 3a 2 (16 + a 2 )) a? 2 

+ (16 (16 + a 4 ) + a 2 (64 + 3/3 4 ) - 2a 6 ) a 8 ai 4 + (8 (8 + a 2 ) - a 4 ) aioai4 

- (16 + 3/3 2 ) ai2«i4 + a?4> (A-3) 



* Notice that the formula for i?i 6 given in the next page can only be found in the ArXiv version of the 
manuscript. For shortness it has been removed in J. Phys. A. published version of the manuscript. This is 
however the only difference between the two documents. 
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-^16 = ((a2«4 - «6 - 16 (4 + 3/3 2 ))a6 - (8 + O2) 2 )a8 a io + (16 + «2) + 04) «6«io 
+2 (8 + a 2 ) a 8 a? - <*io - 2§ (64 + /? 2 (16 + 2 (20 + 17/3 2 ))) a 4 ai ai2 
-128((-8 + /3 2 (2 + 7/3 2 )) a\ + 2 (16 + /3 2 (-4 + /3 2 (-3 + 2 2 ))) a & )a w a 12 
-4(16 (-8 + a 2 (14 + 3a 2 )) + (4 + 5a 2 ) a 4 )a 4 a 6 ai ai 2 + agcuioa^ 
-(4 (-160 + (-24 + a 2 ) a 2 ) + (-16 + a 2 ) tt 4 )a\aiQtti 2 
+ (8 + a 2 ) (8 (16 (8 + a 4 ) + a 2 (48 + 2a 2 + a 4 )) + (24 + a 2 ) a 6 ) a 8 a w a 12 
+ (16 (80 + /3 2 (32 + (3 2 (11 + 2 ))) + 2a 2 (12 + a 2 ) a 4 - (40 + 3a 2 ) a 6 )a 2 w a 12 
+4q 6 (96 (80 + a 4 ) + a 2 (2368 - a 2 (-192 + a 4 ) + 12a 4 + a 6 )) a? 2 
+12(8 + a 2 ) 2 (12 + a 2 ) a 8 a 2 2 - (8 + a 2 ) (12 + a 2 ) (24 + a 2 ) a w a 2 12 

-2(2 8 /3 2 ("80 + (-8 + a 2 ) a 2 ) a 4 - a\ (a 6 + 4 (4 + /? 2 ) 2 - (-8 + 2 ) /3 4 ) )a 8 a 14 

+ (32 (64 + (3 2 (2 + /3 2 ) (8 + 5a 2 )) + a 4 (-2 9 + 44a 2 2 + 3a 2 (64 + a 4 ))) a 6 ai ai 4 

- (28 (8 + a 2 ) + 3/9 4 ) «6«ioai4 - 8 (8 + a 2 ) (80 + a 2 (22 + a 2 )) a 8 aioai4 

- (8 + u 2 ) ((24 + a 2 ) a 4 + 4a 6 )a 8 a w a u - 16 (/3 2 (8 + 30 2 (4 + f3 2 )) - 32) a 4 a 2 4 
+ (4 (-32 + 2a 4 + a 6 ) - a 2 (4 (8 + a 2 ) + /3 4 )) a 2 ai 4 

+4 (32 (2 9 + /3 2 (4 + 2 ) (88 + 2 (20 + 9/3 2 ))) + 2 (40 + 7a 2 ) c$ a 12 a u 
+16a 4 (-320 + 3/J 2 (-16 + 2 (20 + 70 2 ))) ai 2 ai 4 

-(4 (896 + 2 (528 + /J 2 (112 + 2 ))) - (-192 + (-32 + a 2 ) a 2 ) a 4 )a 6 ai 2 ai 4 

-(((«2 - 8) «g + 3(8 + a 2 ) 2 (16 + a 2 ) a 8 ) - (8 + a 2 ) (72 + 5a 2 ) ai )ai 2 ai 4 

-16 (2 10 + 2 (896 + 2 (416 + /J 2 (108 + 13/3 2 )))) a 2 4 

+ (2432 + 24a 4 + a 2 (688 + 52a 2 + 5a 4 ) - 2a 6 ) a 6 a 2 4 + 3(8 + a 2 ) 2 a 8 a 2 4 

-4 (8 + a 2 ) ai a 2 4 + 2(8 + a 2 ) 2 ((8 + a 2 ) a 8 - "io) "10016 

+ (8 + a 2 ) 2 (-8 (2 s + a 2 (80 + 8a 2 + a 4 )) + (24 + a 2 ) a 6 ) a 12 a w 

+ (8 + a 2 ) 2 (1152 + 8a 4 + a 2 (336 + 28a 2 + 3a 4 ) — 4ag) a^aiQ — (8 + a 2 ) 4 a 2 6 

-2 10 (896 + a 2 (496 + a 2 (100 + 7a 2 ))) aioai 2 /3 2 - 32a 4 1 /? 2 (2a 10 ai2 + a 8 a M ) 

- (16 + 3a 2 ) a\a\ 4 2 + 2 lA a 8 a lA 2 2 (2 + /3 2 ) - 2 s a\a 8 ai A (-8 + a 2 /? 2 ) 
+64 (-2 8 (88 + 3a 4 ) - 2a 2 (7936 + 2a 2 (1120 + a 2 (144 + 7a 2 )))) a\ 2 
+64 ((320 + a 2 (72 + 5a 2 )) a 4 - (6 + a 2 ) a 2 /3 2 ) a\ 2 + 12a 2 /3 2 ai ai 4 /3 4 

+ (-16 (16 + a 2 (9 + a 2 )) a 4 - a\0 2 - 64 (16 + 2 (8 + /3 2 (4 + /? 2 ) ))) a? 
-16 ((2 + a 2 ) a 4 - 8 (28 + 5a 2 )) a 6 a 8 a l4 4 + 4 ((a 4 - 32) /? 2 - 64) 4 a 8 a 2 lo 
+4 (4a 4 (/3 2 (12 + 110 2 ) - 32) + 32 (64 + 2 (32 + 2 (18 + 5/3 2 )))) ai a 14 /? 4 
+8 (-4 (48 + a A ) + a 2 (-32 + 2a 2 + a 4 ) - a 6 ) a§a 8 ai 2 + 2(8 + a 2 ) 2 a 2 8 a u 6 
+64 ((2 (-48 + a 4 ) + a 2 (-8 + 2a 2 + a 4 ))) /3 4 a 6 a 8 ai 2 - 8(8 + a 2 fa 2 8 a l2 & 

- (8 + a 2 ) (8 (48 + a 4 ) + a 2 (16 (9 + a 2 ) + 3a 4 ) - 4a 6 ) ai ai 6 /? 6 

-(8 + a 2 ) 2 a 8 a 16 0l + (-a 6 + 2 4 ) 0% (ia 2 10 + (16c*i 2 - 8a u + «i 6 ) 6 ) . (A.4) 

The polynomials R 2n can be obtained from the polynomials R 2n _ 2 by substituting au — > a^+i, 
/3 fc — > 0k + i, multiplying by a global factor —1 and introducing as many factors o\ as needed 
in order to achieve that every term in R 2n is a product of exactly n — 3 elementary symmetric 
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polynomials. For example: 



K = o, 

Rg = a i, 

R°io = -<*l (A.5) 
R°2 = — (8(2«i + 03) — 05)0:107 — [5 3 {a\ — 2oiOg) + a\ on, 

R° 4 = f3 3 a 7 ((3 5 a 7 - ((16ai + o 3 )a 9 - 3aian)) + aiaf + /3§(8(2/? 3 + o 5 )o 9 - 2/3 5 an + P 3 a 13 ), 
and so on. Notice that 01 = o\. 
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